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ABSTRACT 

Using a novel technique, we achieve ^100 picoarcsecond resolution and set an upper bound of less 
than 4 km for the characteristic size of the Vela pulsar's emission region. Specifically, we analyze 
flux-density statistics of the Vela pulsar at 760 MHz. Because the pulsar exhibits strong diffractive 
scintillation, these statistics convey information about the spatial extent of the radio emission region. 
We measure both a characteristic size of the emission region and the emission sizes for individual 
pulses. Our results imply that the radio emission altitude for the Vela pulsar at this frequency is less 
than 340 km. 

Keywords: Methods: data analysis - Methods: statistical - pulsars: individual Vela pulsar - Scat- 
tering - Techniques: high angular resolution 
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1. INTRODUCTION 

1.1. Pulsar Emission and Scintillation 

Owing to their extraordinary emission, pulsars are in- 
valuable probes of a vast array of extreme physics, from 
the supranuclear density interiors and the ~10 12 G mag- 
netospheres of neutron stars, to the turbulent plasma 
of the interstellar medium (ISM). The regularity of 
their time- averaged emission has been applied to sensi- 
tive tests of general relativity (jTavlor fc Weisberall989t 
IKramer et a.1.1 [2006h . and the achieved timing precision 
holds promise for direct observation of gravitational 
wave s (|Detweilerl Il979t Uenet et all 120061: iHobbs et all 
120101 ). Yet, remarkably, no consensus exists for the 
precise origin of this emission; its geometry, stabil- 
ity, and physical mechanism remain e nigmatic , even 
after nearly a hal f-century of study (Melrose l2004t 
iVerbiest et al.l[2~009f ). The incredible compactness of the 
pulsar and its surrounding plasma prohibit traditional 
imaging of the emission process, but a number of meth- 
ods can achieve the requisite resolutio n, using the in- 
terstellar pla s ma as an effective ' lens' (|Lovelacej 119701: 
Backer! 119751: ICordes et al.l 119831: iCornwell fc Naravanl 
19931: iGwinn et alJll998t lShishovll2oTo[ )~ 



Variations of this stochastic lens give rise to the ob- 
served scintillation, and statistics of these variations con- 
vey information about the spatial structure of the emis- 
sion. Longer observing wavelengths enhance the signa- 
ture of this structure because the stronger scattering ef- 
fectively enlarges the aperture. However, the increased 
spatial resolution comes at the expense of spectral reso- 
lution, and it becomes difficult to decouple the scintilla- 
tion from source noise and variability. Indeed, the pulse- 
broadening timescale can easily exceed that of intrinsic 
pulsar variation. 
We can circumvent these difficulties by constructing 



spectra over accumulation times that include all pulsed 
power. Intrinsic sub-pulse variability then affects the cor- 
relation of spectral noise but not the mean spect rum or 
single-channel statistics (jGwinn fc Johnsoril20llT ). Fur- 
thermore, we can distinguish the effects of extended 
source geometry from the scattering evolution and pulsar 
variability by varying the degree of temporal averaging; 
with no averaging, only the source emission structure 
contributes. Additionally, we can identify the emission 
sizes of individual pulses. We have p resented a detailed 
math ematical treatment separately ([Johnson fe Gwirml 
l20ll hereafter JG12). 

We constructed spectra in this way and then com- 
pared histograms formed from the spectral samples with 
appropriate models for a strongly scintillating source. 
Such models benefit from minimal assumptions about the 
nature and distribution of both the scattering material 
and the underlying emission; this generality is achieved 
at the expense of resolution in pulse phase. The in- 
ferred emission size is therefore a characteristic size of the 
emission region, quantified via a dimensionless parame- 
ter 7s, which depends on the spatial standard deviation 
of source intensity, weighted by integrated flux density. 
Translating 7 S to physical units at the pulsar additionally 
requires specification of the scattering geometry. 

1.2. Comparison with Previous Work 

Several previous investigations have also quantified 
scintillation statistics in order to estimate the size of the 
Vela pulsar's emission region. In general, for a strongly 
scattered point source, the scintillation acts upon in- 
tensity as a multiplicative gain Q that is the squared 
modulus of a complex random walk. The probability 
density function (PDF) of Q is therefore exponential, 



P{Q) 



e , with an associated modulation index of 
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unity: m 2 = {G 2 ) / (G) —1 = 1. A spatially-extended 
emission region smoothes this variation, and so the mod- 
ified P{G) or reduced modulation index can be used to 
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estimate the emission size (jSalpeterl [l967t iCohen et al.l 
[19671 IGwinn et al.lH998l) . 

The principal difficulty in detecting these modifications 
is the additional noise that arises from the noiselike na- 
ture of the source, known as self-noise, and from the 
background. Pulsars are particularly challenging targets 
because their dramatic intrapulse and interpulse varia- 
tions also contribute modulation. All of these sources of 
noise can easily imitate or mask the signature of an ex- 
tended emission region. Previous efforts have attempted 
to minimize these sources of noise by spectral averag- 
ing in both frequency and time. In contrast, the present 
technique fully accounts for the effects of self-noise and 
intrinsic pulsar variability. 

For example, IGwinn et al.l (|1997l) sought to detect the 
modification of P(G) through the PDF of the correla- 
tion function amplitude on the Tidbinbilla-Parkes base- 
line at 2.3 GHz. Because this baseline was much shorter 
than the scale of the diffraction pattern, the observa- 
tional setup was effectively a zero-baseline interferome- 
ter. Their best-fitting model corresponde d to an emis- 
sion re gion with a FWHM of 460 ±110 km. IGwinn et all 
(2000) then extended this analysis to three gates across 
the pulse, also at 2.3 GHz. They improved the model 
PDF and quantified the effects of their averaging in fre- 
quency and time. They estimated a declining size of the 
emission region from 440 ± 90 km to less than 200 km 
across the pulse. 

Rather than anal yzing distribution functions, 
iMacquart et al.l ([2000) used the modulation index 
to quantify the effects of emission size for single-dish 
data. They measured m 2 w 0.9 at 660 MHz and thereby 
inferred the extent of source emission to be no more than 
50 km (FWHM = 120 km). Including the contribution 
of self-noise leads to an in ferred FWHM of the emission 
region of -200 km OGTll. 

More recently, IGwinn et al.l (|2012l ) analyzed visibilities 
on the Tidbinbilla-Mopra baseline at 1.65 GHz and fit 
models to the projections of the real part and imaginary 
variance. They found the size of the emission region 
to be large both early and late in the pulse (—400 km 
and —800 km, respectively) but near zero in the central 
portion, where the pulse is strong. 

An analysis at lower frequencies has the advantage of 
increasing the signature of size but the disadvantage of 
incurring more mixing between the scintillation and the 
pulsar noise. We demonstrate that this difficulty is not 
prohibitive, and that a careful statistical description of 
observed power spectra at low frequencies can yield fresh 
insight into the physics of both interstellar scattering and 
pulsar emission. 

1.3. Outline of Paper 

In 521 we briefly review the theoretical descriptions 
of pulsar emission and interstellar scintillation. Then, 
in 53 we describe our observations and the subsequent 
construction of both the observed and model PDFs of in- 
tensity; we also present the expected residual structure 
arising from an extended emission region and outline the 
model fitting procedure. Next, in 531 we give our analysis 
for various pulse selection cuts and derive our estimates 
for the size of the emission region. In 53 we quantify 
possible biases and systematic errors in our inferences of 
the size of the emission region. Then, in 53 we interpret 



our measurements in terms of beamed, dipolar emission, 
and derive the corresponding emission altitudes. Finally, 
in 53 we summarize our results, their impact on our un- 
derstanding of pulsar magnetospheres, and the prospects 
for future work. 

2. THEORETICAL BACKGROUND 

2.1. Pulsar Radio Emission 

A full understanding of pulsar emission is elu- 
sive, but proposed scenarios typically invoke a rel- 
ativistic electron-positron plasma that flows outward 
from the pulsar polar caps to the light cylinder 
(jGoldreich fe Julian]|1969ft . This outflow then generates 
radio emission via coherent curvature r adiation (|Sturrockl 
[1971 : IRuderman fc Sutherland! fl975h . The emission 
is beamed, ducted, and refracted as it propagates 
through the upper magnetosphere dBarnard fe Aronsl 
1986; Arons fe Barnardlfl98i ILvutikov fe Parikhl l2000'). 

The precise altitude at which the emission is 
generated is uncertain, and various analyses sug- 
gest emission heights fro m a few stellar radii 
( Arons fc Scharlem ann 1979) to hundreds of stellar radii 
(Karastergiou fc Johnstod[200l for Vela-like pulsars. In 
the classification and interpretation of lRankinl (|1990l ) , the 
Vela pulsar is a core-single (St) pulsar, with emission 
arising effectively at the stellar surface. 

2.2. Interstellar Scattering 
2.2.1. Scintillation Physics 

Density inhomogeneities in the dilute, turbulent 
plasma of the ISM scatter the radio emission, leading 
to multipath propagation. At 760 MHz, the Vela pul- 
sar is deep within the regime of 'strong' scattering, for 
which the Fresnel scale is much larger than the diffrac- 
tive scale. Furthermore, the timescale for evolution of 
the diffraction pattern at the observer is several sec- 
onds, so each scintillation element is effectively static for 
the duration of a single pulse (—5 ms). The scattering 
then acts to convolve the pulsar signal with a stochastic 
propagation kernel, w hich reflects the scintillation and 
temporal broadening (lHankinsI 119711 : iWilliamsonl 119721 : 
Gwinn fc Joh nson 201l|). 

The spatial extent of the emission region influences 
the scintillation. A pair of dimensionless parameters, 
{7s. i,7s, 2}, quantify the emission region and give its 
squared characteristic size in an orthogonal, transverse 
basis jx.i , X2} , in u nits of the magnified diffractive scale 
(jGwinn et al.lfl99l : 

7M = ( Jj^TY ) • 

\ 1 271- 0i J 

Here, D is the characteristic observer-scatterer distance, 
R is the characteristic source-scatterer distance, A is the 
observing wavelength, 6*.; is the angular size of the scat- 
tering disk along Xi, and oi is the standard deviation of 
the distribution of source intensity along Xi. This rep- 
resentation accommodates anisotropy of both the source 
emission and the distribution of scattering material. 

Nevertheless, the effects of anisotropy in the emission 
and scattering only manifest at quadratic order in 7 S ; 
thus, we use 7 S = (7 Sj i + j s ^)/2 to avoid parameter de- 
generacy. Likewise, we state our results in terms of the 
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standard deviation a c = <j\ = o~2 of a circular Gaus- 
sian intensity profile. The FWHM of such a profile is 
2 V / In4cr c « 2.35cr c . 

2.2.2. Interpretation of the Inferred Emission Size 

The emission size that we infer from 7 S reflects the 
full extent of the emission region. Explicitly, we observe 
the spatial standard deviation of the emitting region, 
weighted by integrated flux density (|JG12D . Thus, even 
if the instantaneous emission is pointlike, if the received 
radiation arises from spatially offset emission sites then 
7s > 0. We derive the effective emission size of beamed, 
dipolar emission in $6] 

2.2.3. Scattering Geometry of the Vela Pulsar 

The conversion of 7 S to a physical size at the pulsar 
requires measurements of the location, sca le, and nature 
of the scattering. Using VLBI parallax, IDodson et al.l 
(2003) estimated a distance of 287**7 pc to the Vela 
pulsar. Combined with the angular and temporal broad- 
ening of the pulsar, this distance identifie s the charac- 
terist ic location of the sca ttering material (jGwinn et al.l 
f993). iGwinn et al.l (|1997f ) measured the angular broad- 
ening to be (3.3 ± 0.2 mas) x (2.0 ± O.f mas) (FWHM) 
at 2.8 GHz. Combined with the temporal broadening, 
as calculated via the decorrelation bandwidth Av^ = 
66 kHz at 2.3 GHz, these measurements suggest a char- 
acteristic fractional distance of D ^_ R ~ 0.72 to the scat- 
teri ng material and an effective magnification of D/R 
2.9 ()Gwinn et a l. 2000). Thus, if the source emits a cir- 
cular Gaussian intensity profile, then a c ~ ( soo M mhz ) x 
440^/tI km. 

Different models of the scattering medium or emission 
geometry will update this conversion but will not affect 
our model PDFs or measurements of j s . For example, 
we quantify the effect of the phase structure function of 
the scattering material in i)5.3l 

3. OBSERVATION AND DATA REDUCTION 
3.1. Observation and Baseband Recording 

We observed the Vela pulsar for one hour on 22 Oc- 
tober 2010 and one hour on 25 October 2010 using the 
Robert C. Byrd Green Bank Telescope (GBT) and the 
recently built Green Bank Ultimate Pulsar Processing 
Instrument (guppi). For flux calibration, we also ob- 
served 3C 190 (0758+143) for approximately 20 minutes 
prior to each pulsar observation. 

We used guppi in 8-bit, 200 MHz mode, and our ob- 
servation spanned 723.125—923.125 MHz. Rather than 
utilizing the standard real-time dedispersion of GUPPI, 
we recorded the channelized complex baseband voltages 
directly to disk for offline processing. The collected data 
constitute approximately 8 TB. 

3.2. Coherent Dedispersion and Formation of Spectra 

We processed the baseband data using the 
digital signal proce s sing software library DSPSR 
(|van Straten fc Bailesl 1201 1[ ). which first performed 
coherent dedispersion and then formed dynamic spectra 
via a phase-locked filterbank. This technique constructs 
each spectrum from a fixed-length time-series beginning 



at the corresponding pulse phase, rather than window- 
ing the pulse into a fixed number of temporal phase 
bins. Our spectra had 262,144 channels per 25 MHz, 
corresponding to 95.4 Hz resolution. The ^10.5 ms 
spectral accumulation time easily contained all of the 
power in one pulse. 

We exported the coherently-dcdispersed spectra from 
DSPSR into PSRFITS format and then performed flux cal- 
ibration fo r each polarization using the software suite 
psrchive (jHotan et al.l 120041 ). Subsequent processing 
utilized the psrchive libraries, which we integrated into 
our custom software. 

Although reversing the Faraday rotation is straight- 
forward, we chose to analyze the measured linear polar- 
ization streams without modification in order to avoid 
introducing possible artifacts from calibration errors and 
mixing. Consistency of our results between the two po- 
larizations is, therefore, a meaningful and encouraging 
indication that our inferences reflect properties of the 
pulsar rather than instrumental limitations. 

3.3. Data Excision and Bandpass Correction 

3.3.1. RFI Excision 

We iteratively excised RFI using a median smoothed 
difference criterion, similar to that implemented in 
psrchive. Specifically, we compared the estimated back- 
ground noise in each frequency channel with the noise in 
its surrounding 250 kHz sub-band; we rejected channels 
that differed from the sub-band median by more than 
three times the sub-band standard deviation. Such re- 
jections typically constituted a tiny fraction of the data 
(~0.3%, as expected from the 3a cut). Also, to reject 
impulsive RFI, we dropped all pulses for which the off- 
pulse mean differed from a three-minute moving average 
of the off-pulse region by more than 5%. 

3.3.2. Background Noise and Bandpass Shape Estimation 

GUPPI divided the observed 200 MHz frequency band 
into 32 filterbank sub-bands, each of 6.25 MHz. Because 
of spectral roll-off, we only analyzed the central 3.2 MHz 
of each filterbank sub-band; each such segment was ana- 
lyzed independently. 

We estimated both the background noise level and the 
bandpass shape by analyzing the off-pulse spectra. The 
off-pulse region excluded the 10.5 ms on-pulse spectra, 
as well as a frequency-dependent phase range of width 
0.15 that contained negatively-dedispersed power leaked 
from the interleaved sampling. 

Next, we normalized the bandpass. We determined 
the bandpass shape for each three-minute block using a 
(3 minute) x (100 kHz) moving average of the off-pulse 
data, or ~2xl0 6 samples for each spectral channel. The 
bandpass variations across each 3.2 MHz sub-band were 
typically a few percent; thus, they contributed minimally 
to the observed statistics because of their small variance. 

3.3.3. Rejection of Pulses with Quantization Saturation 

Because the Vela pulsar is both bright and variable, 
particularly when coupled with scintillation, some of the 
strongest pulses saturated our quantization thresholds. 
Moreover, because the quantization levels were fixed 
across the band, the strong linear polarization of the pul- 
sar (combined with Faraday rotation) and the receiver 
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bandpass shape caused the prevalence of saturation to 
vary substantially across the 200 MHz observing band- 
width. 

Accordingly, we implemented a five-level 'guard zone' 
at the positive and negative thresholds (of 256 total lev- 
els); we excluded any pulse that included a quantized 
value |x| > 122. This criterion rejected > 75% of pulses 
for the stronger linear polarization near 850 MHz, but 
rejected < 1% of pulses, regardless of polarization, at fre- 
quencies < 775 MHz. To avoid a selection bias against 
strong pulses, we preferentially draw later examples from 
the lower-frequency sub-bands. 

3.4. Observed PDF of Intensity 
3.4.1. Histogram Formation 

We constructed a histogram of measured intensities for 
every (1 hour) x (3.2 MHz) block of on-pulse data (~10 9 
spectral samples). Histogram bins had uniform widths of 
0.05, in units where the mean background intensity was 
unity because of the bandpass normalization. We then 
normalized each histogram by the total number of sam- 
ples and the bin width; we will refer to this normalized 
histogram as the 'observed PDF' of intensity. 

3.4.2. Effects of Finite Histogram Resolution 

In order to compare a theoretical PDF P(I) with an 
observed PDF, we must account for the finite bin width 
of the histogram. Specifically, the theoretical (normal- 
ized) histogram value H (I; AI) for a bin of width AI 
centered on / is given by 



H(I;AI) = 



1 

a7 



I+AI/2 



I-AI/2 



P(l')dl'. 



(2) 



Expanding P(I') about / and replacing the second 
derivative with its discrete representation gives 



AI) « P(I) 



P(I + AI) - 2P(I) + P(I - AI) 
24 



(3) 

This approximation corrects for the finite bin width to 
quadratic order in AI, while only requiring calculation 
of the theoretical PDF at the histogram bin centers. For 
our data, the next order correction, Ai r4 P( 4 )(/)/1920, 
was well below our Poisson noise. The effects of the finite 
bin width were small when N — 1 (typically < 0.005% 
of the PDF amplitude) , but they became appreciable for 
larger N. 

3.5. Model PDF of Intensity 

We independently constructed observed PDFs for each 
scalar electric field (i.e. each linear polarization) and 
with various degrees of temporal averaging N, band- 
width, and pulse selection cuts. In all of these cases, 
the mean source and background intensities of the pulses 
parametrized the point-source model PDF of intensity; 
this parametrization accounted for the intrinsic pulse-to- 
pulse variability as well as the quasistatic evolution of 
receiver gain and noise levels. Because this treatment 
incorporated the observed set of pulse amplitudes rather 
than appealing to assumptions about their distribution, 
our models were robust to aribitrary pulse-to-pulse vari- 
ations, including strongly correlated behavior such as 



nulling or heavy-tailed distributio ns of pulse amplitude , 
as are observed fo r some pulsars ()Lundgren et al.l fl995t 
ICairns et al1l200l . 

3.5.1. Model PDF Parameter Estimation 

We estimated the single-pulse source and background 
amplitudes using the mean of all the on-pulse or off- 
pulse channels in the analyzed sub-band that had not 
been flagged for RFI. For a 3.2 MHz sub-band, this es- 
timate included ^30,000 samples for the on-pulse re- 
gion and ^200,000 samples for the off-pulse region of 
each pulse. This average effectively isolated the intrin- 
sic pulsed power from scintillation because each averaged 
sub-band contained many independent scintillation ele- 
ments (Az/d ~ 1 kHz). However, these estimates were 
subject to both bias and noise, and we determine the 
effects of such errors in <J51 

3.5.2. Construction of the Model PDF 

We now describe the basis and formation of our model 
PDFs of intensity. First, consider N consecutive pulses 
with gated source intensities A\I B , . . . , A^I B and mean 
background intensity /„. After averaging the N s pectra 
formed from these pulses, the PDF of intensity is (| JG12t 
Eq. 3) 

roo —G 

P(I;N)=N dG - (4) 
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The integral over G arises because we assume that each 
spectrum effectively explores the full scintillation ensem- 
ble. 

We calculated this theoretical PDF for every sequence 
of N pulses that contributed to an observed PDF, using 
the estimated model parameters {Ajl s , I n }. We then 
averaged these distributions to yield the 'model PDF' of 
intensity. Note that the model PDF required no fitted 
parameters and was independently constructed for each 
observed PDF. 

3.6. Expected Residual Structure 

The PDF residual 7Z(I) (i.e. the observed minus model 
PDF) contains signatures of extended emission size, 
decorrelation of the scintillation within the averaging 
time, and errors in the model parameters {Ajl s , I n }. 
These perturbations have nearly identical shape, so a 
single fitted parameter quantifies their combined effect. 
We can break this degeneracy by varying the reduction 
parameters such as N and the analyzed bandwidth. For 
example, the temporal decorrelation of the scintillation 
pattern does not contribute when N = 1. 

We approximate the expected cumulative residual aris- 
ing fromthe emission size and the temporal decorrelation 
as (|JGH Eq. 9, 10, 18) 
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This general form expresses the influence of the temporal 
autocorrelation function of the scintillation pattern: 

r _ ([/(ti,/)-/n(ti,/)][/fe,/)-/nfe-,/)]) 

(i(tj)-i n (tj)y 

(6) 

We will quantify alternative sources of residual structure, 
such as parameter errors, in terms of their incurred bias 
on 7 S . As was required for the model PDF, this model 
residual was calculated for every sequence of N pulses 
that contributed to the corresponding observed PDF and 
subsequently averaged. 

If the TV averaged intensities are drawn from consec- 
utive pulses, then the scale TZq of the residual structure 
can be approximated as 



i i i i I i i i — i I i i i i I i — r~n — I i i i i 



AN(N + 1) 
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where Atj 3> N is the decorrelation timescale in pulse 
periods, and we have assumed the standard square-law 
autocorrelation — e~"- % ~i'l ATd l . Thus, the relative 
influence of temporal decorrelation to emission size is 
approximately quadratic in N. 

3.7. Fitting the Residual Structure 

For each observed PDF, we calculated the correspond- 
ing model PDF and residual shape. We then performed a 
weighted least-squares fit to the measured residual after 
correcting for the finite histogram resolution; we assigned 
weights by assuming that the histogram errors reflected 
Poisson noise. The only fitted parameter was the am- 
plitude scale Hq of the residual. In some cases, we per- 
formed multiple reductions on the same set of pulses by 
varying the averaging (N) or the lag between the aver- 
aged pulses (At); we then fit all these reductions simulta- 
neously to the pair of parameters {7 S , Aid} by assuming 
a square-law autocorrelation structure. 

4. ANALYSIS 

We now describe our analysis and subsequent infer- 
ence, in successively tighter pulse cuts. We began by ana- 
lyzing statistics of each polarization, when pulses were all 
analyzed concurrently ( 34.1[) . Then, we analyzed PDFs 
formed by averaging pairs of intensities with varying tem- 
poral lag, to demonstrate that our description of the ef- 
fects of temporal decorrelation is correct ( §4.2)1 . Next, 
we segregated pulses into quantiles by pulse strength to 
investigate the possible covariance between the size of 
the emission region and pulse shape ( ^4.3)1 . Finally, we 
examined the statistics of each individual pulse ( §4.4)) . 

We did not detect the effects of a finite emission region 
in any of these cases, regardless of frequency or polariza- 
tion. Thus, our chosen examples are merely those that 
yield the tightest constraints on 7 S . For this purpose, we 
favor one filterbank sub-band centered on 757.5 MHz, for 
which the polarization position angle of the pulsar was 
nearly aligned with that of the receiver, the pulsar rarely 
saturated the quantization thresholds, the leaked power 
from the interleaved sampling was far off-pulse, and there 
was no evidence for substantial RFI contamination. To 
achieve our tightest limit on the size of the emission re- 
gion, we jointly analyze this sub-band with the adjacent 
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Figure 1. Observed and model PDFs of intensity for three de- 
grees of averaging: N = 1, 3, and 5. There are no fitted pa- 
rameters. The mean intensities of the source and background are 
I S (A) = 0.87 and (7 n ) = 1.00, respectively; the plot therefore ex- 
tends to ~27 times the mean signal. 

sub-bands. The full frequency evolution of our observa- 
tions is a powerful probe of the scattering, on the o ther 
hand, as we analyze separately (| Johnson et al.ll2~012ft . 

4.1. All Pulses 

To utilize the full statistical content of our observa- 
tions, we first analyzed all pulses concurrently, in groups 
only determined by frequency and polarization. We con- 
structed each observed PDF from a (1 hour)x(3.2 MHz) 
on-pulse block of data: ~10 9 spectral samples and 
^40,000 pulses. The scintillation pattern evolved sub- 
stantially over both the span of time and frequency; each 
observed PDF contained samples from ~10 6 independent 
scintillation elements. 

Figure [1] demonstrates one example of the compari- 
son of our model and observed PDFs of intensity; the 
plotted data correspond to the stronger linear polariza- 
tion for the 3.2 MHz sub-band centered on 757.5 MHz. 
The agreement extends over the nearly eight decades 
spanned. The tail of the distribution arises from the 
combined effects of intrinsic variability, noise, and scintil- 
lation; the success of the model indicates that it correctly 
accounts for each effect. 

Our models also correctly characterize the effects of 
temporal averaging, N. The decorrelation timescale in 
pulse periods is At^ ~ 50, so the correction for temporal 
decorrelation is quite small for the plotted reductions: 
7V=1, 3, and 5 (see Eq.E}. 

For N = 1, the residual structure in the stronger polar- 
ization has an amplitude of <10~ 4 , which is comparable 
to the Poisson noise in the histogram. A weighted lin- 
ear least-squares fit to the lowest 200 samples (I < 10) 
in the observed PDF estimates 7 S = (2.0 ± 1.0) xl0~ 4 
with reduced chi-squared x 2 c d = 1-13; the standard de- 
viation of Xred i s a/2/200 = 0.1. Furthermore, because 
intrinsic modulation correlates noise in adjacent spectral 
channels, Poisson weights underestimate the histogram 
errors, so x 2 > 1 is expected. Thus, a conservative upper- 
bound on the effects of size from this single sub-band is 
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Figure 2. Observed and model PDFs of intensity with their 
corresponding residual structure; N = 2. Three cases are plotted, 
corresponding to averaged pairs of intensities that are separated by 
At = 1, 10, and 20 pulses. The residuals are concurrently fit to the 
decorrelation timescale At^ and dimensionless size 7 S . The inset 
shows the fitted residual scale TZo with 3a errorbars as a function 
of At, along with the best-fitting prediction (see Eqs. |5]and[7]>. 

7s < 5 x 10~ 4 , which is equivalent to <r c ss 8 km at the 
pulsar with the conversion of §2.2.31 

For the adjacent sub-bands (centered on 751.25 MHz 
and 763.75 MHz), the fits give 7s = (-2.2 ± l.l)xl(T 4 
and 7 S = (— 1.2±1.0)xlCP 4 . Combining these results 
refines the 3cr size upper-bound to approximately 4 km 
at the pulsar. Narrower sub-bands give similar results, 
with larger errors because of the increased Poisson noise 
and a larger bias from the finite bandwidth (see £)5.1[) . 
In all cases, our results suggest a characteristic size of 
<10 km for the emission region. 

4.2. N — 2 with Temporal Lag At 

We next examined effects dominated by temporal 
decorrelation of the scintillation pattern. We constructed 
PDFs with N = 2 that averaged pairs of pulses separated 
by Ar pulse periods; the amplitude of the residual struc- 
ture (Eq. [5]) then depends roughly quadratically on Ar. 
We approximated the autocorrelation function as Gaus- 
sian, r(Ar) = e~( Ar / ATd ) , to simultaneously fit all lags 
to the pair of parameters, {7 S , Ard}. 

The residuals match their expectations quite well; Fig- 
ure [2] shows one example, again from the stronger linear 
polarization for the 3.2 MHz sub-band centered on 757.5 
MHz. This excellent agreement demonstrates the qual- 
ity of our analytical representation of decorrelation of the 
scintillation pattern within the averaging time. The fits 
estimate 7s = (9.3±5.4)xl0~ 5 giving a 3<7 upper-bound 
on the size of the emission region as ~7 km. 

4.3. Quantiles 

Individual pulse profiles for the Vela pul- 
sar show marked dep e ndenc e on their strength 
(jKrishnamohan fc Downs! 11983ft . To examine co- 
variance between pulse strength and emission size, we 
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Figure 3. Observed and model PDFs of intensity for the subsets 
of pulses in the top and bottom decile by pulse intensity. The- 
oretical residual curves are plotted for a point source (the slight 
variation is from finite histogram resolution) and for a source that 
extends over 5% of the magnified diffractive scale (i.e. tr c S3 20 km); 
the latter is clearly inconsistent with the observed statistics. 

segregated the pulses into quantiles by their mean flux 
density and then analyzed each subset independently. 

This approach has several utilities. Because the re- 
solving power improves with the SNR, a subset of strong 
pulses can maximize the potential resolution. Further- 
more, many instrumental biases vary differently with 
pulse strength than their effects on the inferred size do 
(see M5.2[) ; the variation of 7 S between quantile subsets 
can therefore suggest the presence of such a bias. 

We established deciles by pulse strength. Figure [3] 
shows the PDFs and their residual structure for the 
top and bottom decile; the mean gated signal-to-noise 
(Ajl s /I n ) in these deciles is 1.39 and 0.48, respectively. 
The two residuals are comparable, and both suggest an 
emission size that is substantially smaller than the mag- 
nified diffractive scale. The limit is most stringent for 
the pulses in the top decile because of their higher SNR: 
7s = (1.1 ± 1.8) x 10 -4 , or an 11 km upper-bound on the 
characteristic emission size at a 3cr confidence level. For 
the bottom decile, we obtain 7 S = (— 7.9 ± 5.5) x 10~ 4 , 
which gives a 12 km upper-bound at the same confidence. 

4.4. Single-Pulse Statistics 

Lastly, we examined the statistics for individual pulses. 
For each pulse, our model requires only two parameters 
to describe the spectral statistics: the mean intensities 
of the background (7 n ) and gated signal (AI S ). Poisson 
noise generally limits each estimate of emission size to 
a sizeable fraction of the magnified diffractive scale; the 
sta ndard error c57 s in the fitted 7 S is well-approximated 
by (IJGTl Eq. 26) 



^7s 



1.3 



exp 



0.4 + 0.18 (In S) 2 +0.01 (In S)' 



(8) 



Here, JVtot is the total number of points used to estimate 
the observed PDF, and S = Ajl s /I n is the gated signal- 
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Figure 4. Inferred emission sizes of individual pulses as a func- 
tion of their gated signal-to-noise S; both linear polarizations are 
plotted. Because S determines the standard error for each measure- 
ment to excellent accuracy, we omit errorbars and instead show the 
expected ±3<r errors about 7 S = 0. We do not obtain a statistically 
significant detection of emission size for any pulse. 

to- noise ratio . Thus, for a strong pulse (S ~ 2), 7 S can 
be determined to within ±0.006; whereas, for a weak 
pulse (S = 0.1), 7s can only be established to within 
±0.2. Interestingly, for observations of highly variable 
pulsars, the resolution afforded by strong single pulses 
may exceed that of the entire pulse ensemble analyzed 
jointly. 

Figure 0] shows the fitted sizes for ~30,000 individual 
pulses as a function of their gated signal-to-noise; values 
from both linear polarizations are plotted. This sample 
of pulses included several 'micro-giant' pulses, as well as 
ten s of pulses with t he lat e 'bump' component identified 
by Uohnston et al.l (|2001[ ) . We measured no confident 
detection of emission size effects; the inferred values of 
7 S are consistent with the expected standard errors for 
observations of a source with point-like emission, as given 
by Eq. [8] 

5. SOURCES OF BIAS IN THE ESTIMATED 
EMISSION SIZE 

Many practic al limitations can bias the inferred emis- 
sion size (|JG12t §4.4). For instance, the model parame- 
ters Aj and J n are random variates estimated via finite 
averages of on-pulse and off-pulse spectra; they are there- 
fore subject to noise, which is unbiased, and instrumen- 
tal limitations such as quantization, which can introduce 
bias. Both of these types of errors can mimic the effects 
of an extended emission region. Here, we estimate the 
bias to 7s resulting from each of these effects. We also 
account for uncertainties in the detailed distribution and 
nature of the scattering material. 

5.1. Effects of Noise in Model Parameters 

Noise in the model parameters introduces spurious 
modulation. Subsequent fits to the spectral statistics will 
then exhibit a positive bias in 7 S to quench this excess 
modulation. Fortunately, we can readily characterize the 
expected noise SAj and 5I n and can therefore anticipate 
the influence of such a bias. 



For example, the noise SAj in the source amplitudes 
is dominated by the finite sampling of the scintillation 
pattern for a single sp ectrum ; the corresponding bias is 
7s -> 7s + V(4iV 8cint ) (|JGT2l : §4.4.2), where A scint is the 
number of averaged scintillation elements used to esti- 
mate Aj. 

On the other hand, because the background noise is 
white and approximately stationary, SI n only depends on 
the number of samples Ny, that are averaged to estimate 
I n . However, the incurred bias to 7 S also depends on the 
gated signal-to-noise ratio S: 7 S — > 7 S + l/(4S 2 iVb). 

For each 3.2 MHz sub-band, we had iV sc i n t ~ 
(3.2 MHz)/Ai/ d « 3,000, and N h ~ 10 5 . Both biases to 
7s were therefore within our quoted confidence. Further- 
more, each bias scales as the reciprocal of the analyzed 
bandwidth. We observed no evolution in the inferred val- 
ues of 7s for 1.6 MHz and 0.8 MHz sub-bands and thus 
confidently established that this class of parameters error 
did not substantially influence our results. 

5.2. Effects of Bias in Model Parameters 

Instrumental limitations can also lead to errors in the 
model parameters. For example, interleaved sampling re- 
sults in a negatively-dedispersed pulse 'echo'; for highly 
dispersed pulsars, such as the Vela pulsar, this process 
can therefore leak pulsed power into the off-pulse region, 
leading to a positive bias in the estimated background 
noise: I n — > I n + SI n . Alternat ively, quantization intro- 
duces signal -dependent noise (|Jenet fc Anderson! 119981 : 
IGwinnl [2006T ): the off-pulse region then underestimates 
the effective on-pulse noise. 

Of course, any bias in the estimated background noise 
incurs a corresponding bias in the estimated source am- 
plitude, and vice versa: 5I n = —SAjI s . These errors then 
affect the inferred emission size as 7 S — > 7 S — 5I n /(2AI s ). 
Because the dependence of emission size with pulse 
strength may vary differently than the noise bias, seg- 
regating pulses by their gated signal-to-noise can iden- 
tify the influence of this type of bias. Because we did 
not observe variation in 7 S between decile subsets, this 
type of parameter error is unlikely to have significantly 
corrupted our emission size inference. 

We further tested for the presence of a noise bias by 
comparing the estimates of I n as a function of pulse 
strength. Leaked power, for instance, would have in- 
troduced a bias that was positively correlated with pulse 
strength. We detected a slight negative trend in our es- 
timates of I n with pulse strength; the variation in (7 n ) 
between the top and bottom deciles (as shown in Fig- 
ure [3]) is approximately 0.15%. If this variation reflects 
a bias in the estimated noise, then the incurred contri- 
bution to 7 S is at approximately the 3cr errors we quote. 
However, the covariance between the fitted values of 7 S 
and (/„) is marginal, suggesting that this variation in 
(I n ) is not responsible for any appreciable bias in the 
emission size inference. The variation may instead re- 
flect noise in / n , which weakly anti-correlates with the 
gated signal-to-noise of the corresponding pulse. Then, 
for the full pulse ensemble, the error estimates of £15. H are 
still appropriate, although the segregation into quantiles 
will somewhat exaggerate this effect. 

Overall, these tests suggest that our 8-bit quantization 
scheme and high-quality backend (guppi) effectively con- 
trolled any systematic errors in the model parameters. 
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Remarkably, the determination of 7 S is then limited by 
the Poisson noise in the observed PDF of intensity (con- 
structed from ^10 9 samples). Because the transverse size 
at the source is proportional to ^/t7, our resolving power 
scales as the fourth root of the observing duration and 
bandwidth. 

5.3. Effects of Modified Scattering Assumptions 

Thus far, we have assumed thin-screen scattering 
with a square-law phase structure functio n, as is ob- 
serva ti onally motivated f or the Vela pulsar ([Williamson! 
1972t iKomesaroff et all Il972t iRoberts fc, Ablesl Il982t 
Johnston et al.l 119981 : Uohnson et al.l I2012I ). However , 
even with relaxed scattering assumptions, the same form 
of Eq. [T] still applies, with modified scaling. 

For example, the translation to a general power- 
law structure funct ion -D</,(x) oc |x|" is given by 
(jJohnson et al.ll2012D 



Hence, if the phase structure function is Kolmogorov 
(a = 5/3), then the effects of a small emission size 7 S 1 
will be somewhat more pronounced than for square-law 
scattering. The inferred transverse sizes will therefore 
be smaller if we assume a = 5/3 than those if we as- 
sume a = 2. Our given limits 7s, Q =2 Si 10~ 4 translate 
to 7s,q=5/3 Ss 2.3 x 10 -5 , so this alternative scattering 
environment roughly halves our inferred limits. 

The translation from 7 S to a dimensionful emission size 
also depends on the bulk scattering geometry (see Eq.[T]). 
In this case, the uncertainty in the magnification D/R 
directly affects the inferred emission size. However, the 
excellent determination of both the temporal and angular 
broadening of the Vela pulsar lead to relatively sma ll 
uncertainties for the magnification (jGwinn etaHHooa), 
and thus the uncertainties for the inferred emission size 
are within the standard errors from Poisson noise. 

5.4. Effects of Proper Motion 

Because 7 S represents a characteristic size of the emis- 
sion region, effects such as the proper motion of the pul- 
sar can be important. However, the Vela pulsar has a 
transv erse velocity of only ~60 km/s in its local environ- 
ment (Do dson et al.l [20031 ) - quite slow relative to other 
you ng pulsars and the typ ical neutron star birth veloc- 
ity (|Lvne fe Lorimerlll994f) . Over our 10 ms accumula- 
tion time, the pulsar travels only 600 m. Moreover, since 
size effects arc weighted by pulse amplitude, the relevant 
timescale is the mean width of individual pulses, typically 
~2 ms, contributing an effective displacement of ^100 m 
- much smaller than our present detection threshold. 

6. INTERPRETATION 

The simple, conventional picture of pulsar radio emis- 
sion beamed along lines of a dipolar field provides contact 
between our measurements and those of traditional ap- 
proaches to estimate the emission altitude r om . Specifi- 
cally, this picture enables a translation between the effec- 
tive transverse size of the emission region, er c , and r cm . 
Our measurements of a c imply tight limits on r cm , in 



the context of this simple model, while obviating many 
assumptions of complementary techniques. 

In £16 - H we review the basic geometric model for pulsar 
emission. Then, in £)6.21 we demonstrate how this model 
relates the transverse size of the emission region to the 
emission altitude. Next, in ^6.31 we compare this tech- 
nique and its result with traditional methods and their 
conclusions. Finally, in Q6A[ we discuss possibilities for 
the discrepancy between our result and previous scintil- 
lation analyses of the Vela pulsar. 

6.1. Geometric Model for Pulsar Emission 

A highly-simplified, conventional model describes the 
geometry of pulsar emis s ion at centimeter and longer 
wavelengths. IDvks et al.l (|2004D provide an overview of 
this model, with a discussion of its assumptions and pa- 
rameters. We now outline the salient features for com- 
pleteness. The model assumptions include a dipolar mag- 
netic field, radio-wave emission tangent to field lines, 
and a single emission altitude for identifiable features 
of the average-pulse profile. This model follows from 
theoretical predictions of an ultra-relativistic electron- 
positron wind, streaming alon g the open dipolar field 
lines in the co-rotating frame dGoldreich fe Julian! 1 19691 : 
iRuderman fe Sutherlandlll975l: lMelrosell2004ft . so we as- 
sume that the emission is beamed tangent to the field 
lines into an infinitesimal solid angle. In this case, an 
observer receives radiation from a single point in the pul- 
sar's magnetosphere at each pulse phase. The variations 
of the spatial displacement of this point with pulse phase 
then determine the apparent size of the emission region. 

Because we merely seek a rough correspondence be- 
tween the size and altitude of the emission region, 
we ignore effects of relativistic terms, modifications to 
the dipolar field struc ture, and magnetosphe r ic duct- 
ing and reprocessing (f Blaskie wicz et al.l 119911: iPhillipsI 
1992tlKapoor fe Shukrel|1998HHibschman fe Aronsll2001l: 
Arons fe Barnardlll986l : lLai et alll200lD . 

A pair of parameters define the orientation of the pul- 
sar: the inclination of the observer relative to the rota- 
tional axis (£) and the impact angle of the magnetic pole 
with the line of sight (f3). We also employ standard an- 
gular coordinates: the longitude or rotational phase of 
the pulsar (0) and the angle between the point of emis- 
sion and the magnetic pole (#em (</>))■ For a fixed emis- 
sion altitude r em , the spatial displacement of the point of 
emission -D cm (</>) then follows from the dipolar geometry: 
|-Dem| ~ f em tan _1 (# em /2). The cosine formula of spheri- 
cal trigonometry relates 9 em to pulse longitude <\> and the 
orientation parameters {/3, C}, and the well-known for- 
mula for the swing of position angle with 4> gives the di- 
rection of D (jKome sarofi 19 70l). ILvne fc Graham- Smith 
(|1998t Figure 16.141 and lLorimer fc Kramer! (|2004l; Fis> 
ure 3.4) provide useful figures and discussion. 

The Vela pulsar benefits from several features that per- 
mit confident inferences of the orientation parameters. 
Fits to the X-r ay emission from its pulsa r wind nebula 
yield C » 63° (|Ng fc Romanill200i 120081. In addition, 
the swing of the position angle of polarization across the 
pulse indicates that /3 » —6 .5° ((Radhakri shnan fc Cookd 
Il969t I Johnston etldl 12005ft . Thus, the Vela pulsar is a 
nearly-orthogonal rotator, with the line of sight close to 
the plane of rotation. Note that the duty cycle of the 
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Vela pulsar is quite small. The full angular ext ent of the 
emission is only p w 16° (pohnston et al J 1200 1). and the 
effective angular extent, expressed as the standard de- 
viation in pulse phase weighted by flux density, is only 
p cS = 3.0°, from our measurements. 

6.2. Emission Altitude from Effective Transverse Size 

We now relate the effective emission size o~ c to r cm , 
within the context of this geometrical model. A per- 
fectly orthogonal rotator (f3 = 0) provides an enlighten- 
ing demonstration. In this case, the standard deviation 
of the transverse spatial displacement, weighted by flux- 
density, is <7 c ff ~ Peff^em/3. The factor of 1/3 reflects 
that the angular difference between the tangent and ra- 
dial directions of the dipolar field is roughly one third 
the colatitude of the tangent direction. The translation 
to an effective circular emission region o~ c then incurs an 
additional factor: a c g = V2o~ c . Applying our most strin- 
gent limits for the size of the emission region, o~ c < 4 km, 
we obtain r cm < 3-v/2cr c / p c g fts 320 km. 

Applying the measured geometry for a nearly- 
orthogonal rotator, appropriately weighted by the indi- 
vidual pulse profiles for the stronger linear polarization 
at 757.5 MHz, we find r om w 86er c < 340 km. More- 
over, if we assume that the phase structure function of 
the scattering medium is Kolmogorov (see !j5.3|) . then we 
obtain r cm < 170 km. 



6.3. Comparison with Previous Emission Altitude 
Estimates 

Our results are comparable with the typical conclu- 
sions derived fro m pulse profil e analy sis and polarime- 
try. For example, iKiiak fc Gill (|1997| ) calculated pulsar 
emission heights by measuring pulse profile widths at 0.4 
and 1.4 GHz and assuming that the emission spanned the 
open field-line region at a single altitude. For the Vela 
pulsar at 757.5 MHz, their model predicts r cm w 360 km. 

More generally, the emission may arise from field 
lines that span only a portion of the open field-line re- 
gion. Researchers typically characterize this possibil- 
ity with the fractional colatitude of active field lines 
fg ?s #max y Rlc /fcm- The inferred emission altitude is 
quite sensitive to the assumed emission structure, scal- 
ing as ^fg 2 - However, with strict symmetry assump- 
tions, fg can be inferred in p ulsars with clearly delineated 
core and cone co mponents (|Gupta fe Gangadharall200"3l : 
iDvks et al"1l2004D . Such measurements suggest that fg 
generally lies within the range 0.2—0.8 for conal emis- 
sion. 

Alternatively, the aberration al shift in position angle 
indicates the emission height dBlaskiewicz etall 119911 : 
iPhillipsI [19921: iHibschman fc Aronsl 120011) . This shift is 
insensitive to fg and thus gives a robust r cm without 
constraining fg. For the Vela pulsar, the position a ngle 
offset indicates r cm ~ 100 km ([Johnston et al.lfeOOll ). 

A measurement of o~ c concurrently determines fg and 
r cm , regardless of the profile morphology or the scatter- 
broadening. Thus, it can probe the emission structure 
at low frequencies, where both the size of the emitting 
region and the effects on the scintillation are most pro- 
nounced. 



6.4. Comparison with 1.66 GHz and 2.3 GHz 
Scintillation Studies 

The sizes that we infer for the emission region are 
substantially smaller than those found by the previous 
scintillation analyses discussed in m.'2\ however, our lack 
of resolution in pulse phase precludes a direct compar- 
ison. For example, if the emitting regions are largest 
when the pulsar is weak, then the effects of emission size 
will be diminished. Indeed, the phase-resolved analysis 
at 1.66 GHz suggests precisely this scenario. Further- 
more, our observations span lower frequencies than pre- 
vious analyses. Although the radius-to-frequency map- 
ping (jRuderman fc Sutherland! 119751: ICordeslll978D sug- 
gests that the emission height decreases with increasing 
frequency, the Vela pulsar possibly does not conform to 
this picture. Alternatively, the pulse profile shows con- 
siderable intrinsic evolution at 1 GHz, and the disparity 
between the measured emission size may indicate the in- 
troduction of a new, offset emission component at these 
higher frequencies. Pending observations should resolve 
this intriguing discrepancy. 

7. SUMMARY 

We have analyzed flux-density statistics of the Vela 
pulsar at 760 MHz. In particular, we examined the ob- 
served PDF of intensity in various cuts determined by 
frequency, polarization, and pulse strength. Based on 
single pulse estimates of the background noise and gated 
signal-to-noise, we constructed model PDFs in order to 
identify the signature arising from a spatially-extended 
emission region. 

The predicted form for point-source emission showed 
excellent agreement with the data (x 2 e d~^)' g°°d to 
~0.01% over the low to moderate intensity regions of 
the PDF, in which the effects of an extended emission 
region are the most pronounced. This type of compari- 
son requires no fitted parameters and arises from minimal 
assumptions about the characteristics of both the scat- 
tering (i.e. strong and slowly evolving) and source emis- 
sion (i.e. amplitude-modulated noise with small spatial 
extent). Our inferences are therefore extremely robust. 

We did not detect any effects arising from an extended 
emission region. We achieved the tightest limit by jointly 
analyzing all pulses in three 3.2 MHz subbands, cen- 
tered on 751.25, 757.5, and 763.75 MHz; this limit corre- 
sponds to a characteristic size of the emission region of 
o~ c < 4 km at the pulsar at 3u confidence. This limit in- 
volved information about the nature and distribution of 
the scattering material. If the phase structure function 
of the scattering material is Kolmogorov, for instance, 
then the limit decreases to er c < 2 km. 

The inferred size, o~ c , reflects the spatial extent of all 
the emitting regions during a pulse and, therefore, also 
gives an upper-bound on the instantaneous emission size, 
which we assume is pointlike in our derivation of the 
emission altitude. Explicitly, er c corresponds to the stan- 
dard deviation of a circular Gaussian distribution of in- 
tensity that yields equivalent scintillation statistics as 
those observed. Conversion to an emission altitude also 
involves the effective angular pulse width, weighted by 
flux density, which provides a natural mapping to lat- 
eral emission structure. We thereby obtained constraints 
on the emission altitude: r cm < 340 km (square-law) or 
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r cm < 170 km (Kolmogorov). 

We also individually tested pulses for the signature of 
an extended emission region, again without a statisti- 
cally significant detection. For the strongest pulses we 
obtained limits < 50 km, and for the majority of pulses 
we obtained limits < 200 km. However, even for the 
strongest pulses, this bound only marginally constrained 
the emission to be generated within the light cylinder. 
Nevertheless, our results demonstrate the application 
and reliability of this technique. 

Our upper bound on the size of the emission region 
improves previous estimates using scintillation statistics 
by over an order of magnitude. Without invoking sym- 
metry constraints or assumptions about the distribution 
of emission within the open field-line region, we achieve 
estimates of the emission altitude that are comparable 
with alternative techniques that analyze profile widths 
and polarization characteristics. Thus, even for incred- 
ibly compact 'core' emission, scintillation statistics now 
have the ability to refine and extend the information con- 
veyed by bulk polarimetry. 

The excellent agreement between our data and mod- 
els demonstrates that our technique can be similarly 
applied to other pulsars without fear of spurious infer- 
ence. In particular, the quality of our data appears to 
have controlled the systematic uncertainties to be within 
the fundamental restriction set by Poisson noise. Thus, 
future observations, perhaps at lower frequencies or of 
other pulsars with favorable scattering geometries, may 
confidently measure a spatially-extended emission region 
and thereby greatly enrich our empirical understanding 
of pulsar radio emission. 
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